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ABSTRACT 

This article presents a study of the evolution of the internal structure and seismic properties expected for low/intermediate mass 
Pre-Main Sequence (PMS) stars. Seismic and non-seismic properties of PMS stars were analysed. This was done using 0.8 to 4.4M0 
stellar models at stages ranging from the end of the Hayashi track up to the Zero- Age Main-Sequence (ZAMS). This research con- 
cludes that, for intermediate-mass stars (M>1.3Mo), diagrams comparing the effective temperature {Teff) against the small separation 
can provide an alternative to Christensen-Dalsgaard (C-D) diagrams. The impact of the metal abundance of intermediate mass stars 
(2.5-4.4Mq) has over their seismic properties is also evaluated. 
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1. Introduction 

The study of the Sun's seismic characteristics allows a better un- 
derstanding of its internal structure (e.g. Antia |2005] l. Yet, for 
other stars it is only possible to observe pulsations associated 
^ with low degree spherical harmonics, constraining the results 
that one can obtain (Bedding & Kjeldsen i2003i ). Still, astero- 
P\| , seismic techniques can be used to infer the internal structure of 
C solar-type stars (e.g. Monteiro et al.. 120021) 

The seismic study of Pre-Main Sequence (PMS) stars (e.g. 
Breger [19721 Marconi & Palla 1998, Ripepi & Marconi 2004) 
,—1 ' has been mainly focused on objects located inside the PMS 
f — Instability Strip (IS). However, in terms of the location and size 
of the outer convective regions, low mass PMS stars, located 
near the Zero Age Main Sequence (ZAMS), resemble solar-type 
. stars. Therefore they may be expected to oscillate. 

Not much work has been done regarding the study of PMS 
H solar-type pulsations (Samadi et al. 2005 Pinheiro et al. 120061 ) . 
. 5t This research goes beyond the study of the seismic properties 
of young stars initiated by Pinheiro et al. (2006) by extending 
Monteiro et al.'s ( |2002| l analysis to PMS stars. By taking into 
account the evolution of the seismic properties of young stars 
we evaluate the possibility of testing PMS evolutionary models 
through the use of solar-type pulsations. 



2. Solar-type pulsations 

The Sun, amongst other solar-type stars, displays a type of oscil- 
lation in which pressure acts as a restoring force. These p-mode 
pulsations, known as solar-type oscillations, are stochastically 
driven by the outer convective layers. The power spectrum of 
solar-type pulsators presents two typical frequency separations, 
known as the large (Av„_/) and small {6v„ i) separation. The first 
one coiTesponds to the difference between frequencies associ- 
ated with oscillations of the same degree (Z) and consecutive 
overtone (n), while the latter corresponds to the difference be- 



tween frequencies v„-\j+2 and v„_/. Tassoul ( 119801 ) showed that 
in the asymptotic regime (n » I) these two separations can be 
written as: 
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where Fj , p and p are the adiabatic exponent, pressure and den- 
sity inside the star. Therefore Av is a measurement of the inverse 
of the time that acoustic waves, with a velocity Cs, take to travel 
from the centre to the stellar surface. According to Kjeldsen & 
Bedding (|19951 l, for a fully ionised star composed of an ideal gas 
with an "average" temperature <T>ocM^/R"*, one gets: 



Av oc -^Jm/R^ cc ^|< p> - 
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i.e. Av depends on a global stellar parameter (<p>). The sound 
speed gradient integral r, (defined on Eq. |2]) is related to the 
asteroseismic ratio r„j-6v„j/ Av^j+i (Roxburgh & Vorontsov 
120031 ). Equation 2 hints that r,- and 6v„ i are sensitive to the sound 
speed gradient (dCs/dr) near the stellar centre. Thus 6v„j is sen- 
sitive to the stellar structure and evolutionary status. Diagrams 
which compare these two frequency separations can be used to 
study solar-type stars (Christensen-Dalsgaard 1984) . These are 
known as Christensen-Dalsgaard (C-D) diagrams. 



3. Evolution of the seismic properties 

Stellar models are required to analyse the evolution of the inter- 
nal structure and seismic properties of PMS stars. PMS models 
produced in preparation for the COROT mission (available 
at http://www.astro.up.pt/corot/models) were used in this study. 



2 



F. J. G. Pinheiro: Seismic evolution of low/intermediate mass PMS stars (RN) 



These models were produced using the CESAM stellar evolu- 
tionary code (Morel 11 9971 1. applying the same procedures and 
assuming the same physical ingredients used in the modelling 
of EK Cep (Marques et al. 2004). These Z=0.02 models have 
masses between 0.8 and 3.2M0. For each mass, 13 models were 
available. These are equally spaced in time, ranging from the 
end of the Hayashi track up to the ZAMS. Figure[T]a shows each 
model's position in the HR diagram. 

The model's seismic parameters (Ay & r,) were estimated 
using the integrals shown in Eqs. 1 and 2. The C-D diagram 
displayed in Fig. [T]b was produced following Roxburgh & 
Vorontsov's approach (2003), i.e. r, was used instead of Sv„j. 
The overlap of different models on the C-D diagram implies that 
this diagram is only useful to analyse stars with less than 1 .3M0. 

3.1. Sound speed gradient integral vs. effective temperature 

By comparing, for each model, the effective temperature (Teff) 
against its r, integral (Fig. |3]a) one breaks the small and large 
separation degeneracy displayed by some of the models (Fig. 
|2]b). This means that unlike C-D diagrams, Te/f vs. r, dia- 
grams can be used to study objects more massive than 1.3M0. 
Yet Monte Carlo simulations point out that, for low mass stars 
(M<1.3M0), Teff vs. r,- diagrams can only achieve the precision 
of C-D diagrams if T^ff is accurately known. For instance, a 
3.5% accuracy in the mass determination is achieved either by 
knowing r, and Av with a 5% precision or by knowing r,- (with 
that accuracy) and Teff with a 50K uncertainty. 

Figure |2] shows that the r,- integral evolves in a similar way 
as the average sound speed gradient computed between 0.25 and 

0.33 stellar radii ^< dCJdr >r/4-r/3= 175^ TT'^')- 
hints at r,'s dependence on the inner sound speed gradient. 



This 



3.2. The large separation and the mean stellar density 

From Eq. 3 we get that Avoc<p>'/^. By relying on the definition 
of Av and <p>'^^(Eqs. 1 & 2) we can write the Av/<p>'^^ ratio 
as: 
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where Kq is an estimate of Ay/< p >'^^ obtained using the 
central density (p(0)) and sound speed (Cs(0)). The normalised 
sound speed (Ic) and density (/p) integrals are independent of 
the model's radius, density and sound speed. 

Figure [3] hints that changes in the density and sound speed 
profiles are partially responsible for the observed variations in 
the Av/<p>'^^ ratio. Such profile changes are due to modifica- 
tions in the relationship between pressure and density resulting 
from changes in the location/size of convective/ionisation re- 
gions that occur during the contraction towards the ZAMS. 

This figure also shows that Av/<p>'^^, Ko, Ic and Ip are, to 
some extent, correlated with the model's effective temperature. 
Indeed, the dispersion of the model's Av/<p>'^^ ratio in Fig.[3]a 
is, at any given temperature, around 2 to 5% Thus knowing a 
star's effective temperature allows one to predict its Av/<p>'''^ 
ratio. 
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Fig. 1. a) PMS evolutionary tracks for 0.8 to 3.2Mq stai's. The 
diamonds correspond to the models described in Sect. 3. b) C-D 
diagram produced using the same models. The arrow indicates 
the direction in which evolution takes place. 



4. The contribution of metallicity 

Intermediate mass PMS stars are known for their 5-Scuti type 
pulsations. Samadi et al. (120051) predicted the amplitudes of 
solar-type pulsations for low mass PMS stars located to the right 
of the IS. The size (in solar radius) of the outer convective re- 
gion of Samadi's PMS model #9 is smaller than that of the 
log(Teff)<3.85 models used in Sect. 3. Therefore, one may be 
tempted to extrapolate Samadi's result for log(Teff)<3.85 stars. 
On the other hand, intermediate mass stars located near the 
ZAMS may not display solar-type pulsations since their outer 
convective layers are smaller (about one order of magnitude) 
than the one of Samadi's PMS model #9. The convective cores of 
A type stars (~2M0) can drive g modes (Browning et al. 120041 & 
Antonello et al. 12006b . However, currently it is unknown whether 
core convection can drive solar-type pulsations like outer con- 
vection does. Nonetheless, such conjectures require an extrapo- 
lation of the study of solar-type pulsations to more massive stars. 

In order to explore the metalicity effect on the seismic 
properties of intermediate mass stars, we selected Z=0.01 and 
Z=0.02 models produced by J. Marques to study the PMS 
5-Scuti star VV Ser (Ripepi et al. 120071 ) . These models have 
masses ranging between 2.5 and 4.4M0 covering, for each mass, 
30 evolutionary stages equally spaced in time. Figure |4] shows 
that the seismic properties of Z=0.02 PMS stars evolve in a 
similar fashion to those of Z=0.01 stars. This happens because 
the internal structure of both types of stars evolves in the same 
way. However quantitatively there are some differences. Near 
the ZAMS, metal-poor stars are denser and have larger sound 
speed gradients than their Z=0.02 counterparts. Therefore their 
r, integral is larger. 
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Fig. 2. a) Evolution of the sound speed gradient integral and 
b) the average sound speed gradient < dCJdr >r/4-r/3 for the 
models displayed at Fig. 1 (diamonds). 



5. Asteroseismic test of PMS models 

The seismic properties of a given star can be used to test stellar 
models that reproduce, within given uncertainties, its position in 
the HR diagram. Figure|5]shows this. In this exercise the target is 
a 2.8M0, Z=0.01 star with log(L/L0)=1.51 and log(T,//)=3.77. 
The models tested here are the same ones used in Sect. 4. In 
this case a 20% uncertainty is assumed in the target's log(L/L0), 
a 1% uncertainty in log(Tfyy) and a 5% accuracy in its seismic 
parameters Av and r,. Figure|5]hints that, as expected from Eq. 3, 
Av puts constraints on the model's density. In the same manner 
we can see how r, puts constraints on the stellar evolutionary 
status. Conversely, models that reproduce the seismic properties 
of a given star can constrain its global stellar parameters (e.g. 

Pinheiro et al. l2()03b . 

According to Baglin et al. ( [200TT l. in a 150 day run COROT 
can achieve a O.lyuHz accuracy in individual frequency deter- 
mination. On the other hand, Av„,/^40//Hz and 5v„ /^6;uHz are 
typical frequency separations of intermediate mass PMS stars 
(Pinheiro et al. 120061 1. This means that COROT can determine 
Av„j and 6v„j with accuracies up to 0.5% and 4%, respectively. 
Table[T]shows that some of the models used above have different 
seismic properties, despite occupying the same position in the 
HR diagram. Since these differences are larger that COROT's 
accuracy, one could use COROT to test these models. 
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Fig. 3. a) Normalization of the large separation (Av) by the 
square root of the mean stellar density <p>'^^, b) Av/<p>'^^ ra- 
tio estimated using the core density and sound speed, c) inverse 
of the profile integral Ic and, d) inverse of the square root of the 
profile integral Ip for the models described in Sect. 3 (diamonds). 



6. Conclusions & future developments 

The seismic characteristics of the PMS models evaluated here fective temperahire. Indeed, in the range of metal abundances 
reflect their internal structure. Therefore, as the models contract analysed here, the knowledge of star's eff^ective temperature 
towards the ZAMS their structure and seismic properties change, (with a lOOK uncertainty) allows us to infer its Av/<p>^'^ ra- 
in general terms, the evolution of Av and r, shows no significant tio with a precision between 2 and 5.5%. 

dependence on the model's mass and metallicity. However these Due to a degeneracy between the seismic properties {6v & 

variations are correlated to some extent with the model's ef- r,) of some models, C-D diagrams are only useful for testing 
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PMS models with less than 1 JM©. This degeneracy is broken by 
taking into account their temperature. Consequently, T^ff vs. r,- 
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Fig. 4. a) Evolution of the Av/< p >'^^ ratio and b) r,- in- 
tegral for diferent stellar models. The black diamonds corre- 
spond to 2.5M0, 3.IM0 and 4.0Mq Z=0.01 models. The grey 
diamonds correspond to 2.5-2. 8M0, 3.1-3.4Mo and 4.0-4.3Mo 
Z=0.02 models. 
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Fig. 5. a) Z-OM and, b) Z=0.02 stellar models that reproduce, 
with 5% accuracy, the large separation Av (black diamonds) and 
the r, integral (white diamonds) of a 2.8M0, Z=0.01 test star 
whose position in the HR diagram is given by the black cross 
(log(Tj.yy)=3.77 & log(L/LQ)=1.51). The grey diamonds corre- 
spond to the models that reproduce both seismic properties with 
a 5% accuracy, while the small squares represent models that do 
not reproduce them. The boxes represent a 1 %x20% uncertainty 
region (in log(Te//) and log(L/L0)) around the black crosses. 
The lines are 2.5 to 4.0Mo PMS evolutionary tracks. 



Table 1. Models that occupy the same position in the HR dia- 
gram and have different seismic characteristics. 



M/Mo 




log(L/Lo) Z age(Myr) 


Ay(^Hz) 


r,(s-') 


2.5 


3.75 


1.16 0.01 1.66 


27.9 


-366.4 


2.5 


3.75 


1.16 0.02 2.46 


28.6 


-454.0 






relative difference 


2.3% 


19.3% 


2.6 


3.86 


1.68 0.01 2.25 


24.1 


-762.4 


2.7 


3.86 


1.68 0.02 2.74 


24.4 


-799.1 






relative difference 


1.3% 


4.6% 


2.6 


4.05 


1.92 0.01 2.82 


53.3 


-1919.7 


3.1 


4.05 


1.92 0.02 2.58 


59.2 


-1662.3 






relative difference 


9.5% 


13.4% 



diagrams can be used to analyse stars more massive than I.3M0. 
On the other hand, the accuracy in individual frequency deter- 
mination that COROT can achieve is enough to test several PMS 
evolutionary models. This result supports the study of young 
solar-type pulsators. 

In the near future we will analyse the effect that stellar pa- 
rameters, such as the mixing and the overshooting, have on the 
seismic properties of low/intermediate mass PMS stars. The seis- 
mic properties expected for each model will be estimated using 
the ADIPLS pulsation code (Christensen-Dals gaard i 1 982i l. 
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